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Abstract. In this talk I concentrate on two topics closely related to the 
understanding of the energy and the system size dependence of the elliptic flow: 
. determination of the elliptic flow (i^) "free" from the effects of non-flow and 

^ ' flow fluctuations, and the role of fluctuations in the initial eccentricity of the 

04 . overlap region. I introduce a new approach for the analysis of the distribution 

in flow vector, dP/dq n , namely, I propose to use the Bessel Transform of this 
. ditribution. I show that the Bessel Transform method is similar to the Lee- 

Yang Zeroes method, and is very transparent in its meaning and applications. 
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■ 1. Introduction 

Anisotropic flow, and, in particular, elliptic flow, plays a very important role in our 
understanding of heavy ion collisions at high energies. Since the discovery of the 
in-plane elliptic flow at AGS [1], and later at CERN SPS [2], the attention to this 
$H ' phenomena has been continuously increasing. The first RHIC measurements [3] by 

the STAR Collaboration made it one of the most important observable at RHIC [4] . 
The large elliptic flow, along with its mass dependence at low transverse momenta [5, 
6] are often used as a strong argument for early thcrmalization of the system; the 
observed [7] constituent quark scaling [4, 8] of elliptic flow strongly suggests that 
the system spends significant time in the deconfined state. 

In order to understand deeper the physics of the elliptic flow and the processes 
governing the evolution of the system, one has to measure and understand elliptic 
flow dependence on the system size and collision energy. Having in mind also 
the dependence on the centrality of the collision makes the problem indeed 'multi- 
dimensional'. In [9] the authors suggested that the elliptic flow follows a simple 
scaling in the initial system eccentricity and the particle density in the transverse 
plane, i^/e oc l/S dN c h/dy (where e is the initial system eccentricity and S is the 
area of the overlap region of two nuclei). If true, all the data taken at different 
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energies, colliding different nuclei, and at different centralities should collapse on 
the same universal curve. Indeed the available data are consistent with such a 
scaling [4, 10]. Unfortunately, the systematic uncertainties in the results are too 
large to conclude on how well the scaling holds. The main difficulties are the 
evaluation/elimination of the so-called non-flow correlations (azimuthal correlations 
not related to the reaction plane orientation) , and effects of flow fluctuations, along 
with uncertainties in the calculation of the initial system eccentricity needed for 
such a plot. 

The non-flow contribution is usually suppressed when one uses multi-particle 
correlations to measure anisotropic flow [11,12]. Unfortunately, the measurements 
of flow via correlations of different number of particles are also different in their 
sensitivity to the event-by-event flow fluctuations [12, 13]. Thus it become difficult to 
disentangle different contributions (non-flow and flow fluctuations) from measured 
azimuthal correlations. It was concluded in [14] that the Lee- Yang Zeroes method 
is the best for suppression of non-flow correlations. Initially I thought that this 
method could be less sensitive to the flow fluctuations. Below I show that this hope 
was over-optimistic, and the results of the Lee- Yang Zeroes method are affected by 
flow fluctuations in a similar way as four-particle cumulant results. 

As it was already mentioned, the flow measurements via correlations of different 
number of particles, such as two- and four-particle correlation results, «2{2} and 
t>2{4}, are sensitive to flow fluctuations in different ways. The contribution of 
fluctuations should be taken out before one tries to explore the scaling v 2 /e oc 
1/S dN ch /dy. Here we discuss a different and more convenient approach, namely to 
rescale the flow results with an appropriate value of eccentricity that already include 
the effect of fluctuations. Such type of eccentricities are studied and presented in 
the second part of this talk. 

2. Fourier and Bessel Transforms of dP/dq x and dP/dq 2 
Consider the distribution in x component of n-th harmonic flow vector 



where M is the multiplicity of particles used to define the flow vector. Let us denote 
the distribution in Q x in the case of no collective flow (v n = 0) to be f n f(x = Q x ). 
In the case of large multiplicities M, this distribution can be approximated by 
Gaussian, but we will not need an explicit form here. In the case of non-zero flow 
the distribution can be written as a superposition of "no-flow" distributions, / n f, 
"shifted" in the direction of flow by appropriate amount (for details, see [16]; the 
necessary condition is \[M 3> 1): 



M 




(1) 
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Now let us consider the Fourier transform of this distribution: 

2^ 



f(k) = (e ik <3*)= I dQ x e lk ^ I d lf nf (Q x - v Mcos(n*)) 



= J^J dQ x e tkQ *f nl (Q x - «Mcos(n*)) 

= / ^e^ Mcos ("*) J dte* kt f nl (t) - J (kvM)f nt (k) (3) 

Remarkably, the flow contribution is completely factorized out. As one expects the 
function f n f(k) to be close to Gaussian (at large M it must approach Gaussian 
according to the Central Limit Theorem) , the zeros of the Fourier transform will be 
determined by the zeroes of the Bessel function J a (kvM), thus providing a method 
of extraction of flow from the distribution in Q x . For example, if fci is the first zero 
of the Fourier transform, 

v = joi/kiM (4) 

where joi ~ 2.045 is the first zero of Jo Bessel function. The above result is the same 
as one would obtain applying the Lee- Yang Zeroes method (using "sum" generating 
function); in fact the relation to the Fourier transform was already pointed out in 
the original paper [14]. 

Even more interesting result one obtains considering two-dimensional distribu- 
tion, d 2 P/dQ x dQy. In this case 

r dp 

= J dQJ (kQ)— ~ J (kvM), (5) 

which, as shown above, is reduced to the Bessel transform of the distribution in 
the magnitude of the flow vector. Note that in this approach (valid in the limit of 
\[M 3> 1) the flow contribution is completely decoupled from all other correlations. 
It happens due to the collective nature of flow. Note also that in the same limit 
one expects the distribution in flow vector to be Gaussian due to the Central Limit 
Theorem, thus explaining why the fitting of the distribution to the form derived 
in [16] (such fits have been used in [1, 12]) is also non-sensitive to non-flow correla- 
tions. Thus in this limit all three methods, the Bessel Transform, Lee- Yang Zeroes, 
and fitting q-distribution become very similar, if not equivalent. 

Fig. 1 (left panel) shows the fit to the q-distribution of the simulated data (four 
million events with 400 particles in each event; 50% events have flow v = 0.04 
and 50% with v = 0.06). Right panel shows the Bessel transform of the same 
distribution. One finds that the results obtained with two methods are indeed very 
close to each other. For a discussion of a small deviation from 0.05, expected for 
this case, see below. 

Though our goal here is to calculate the integrated flow, it is straight-forward 
to generalize the method for the calculation of differential flow as well. For that 
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Fig. 1. (left panel) ^-distribution for simulated events shown together with fit 
function. Dashed line shows the expected distribution in case of v = 0. (right 
panel) Bessel transform of the distribution on the left. (See the last row in Table 1 
for the flow values obtained in this case). 



one considers 

(cos(^)e ifeQ * 



= / ^dcf> i (l + 2vcos(cj> i -*))co S (ci>)e ik Q*f(Q x ) 



2^ 



vcos(yy k6Mcos ^\f nf (k) = v,h{kvM)~f ni {k), 



(6) 



where we have introduced notation v for the average flow in the region of the flow 
vector definition. From the above equation one can obtain the flow v by different 
methods. The one suggested in [14] is based on the relation 



(Q x e ikQ *) = vMJ 1 (kvM)f ni (k). 
Combining Eqns. 6 and 7 one finds: 

v _ g M (coB(0i)e <fc Q-) 

Alternative way would be to use the expression 3 directly. Then 

_ (cos(<fc)e ifc Q-) Jp(kvM) 
v ~ (^fcQT) J^kvM) 



(7) 



(8) 



(9) 



More detailed discussion of the differential flow measurements using Bessel Trans- 
form is deferred for later publication. 

The main advantage of using the Bessel Transform (as well as Lee- Yang Zeroes 
method, or fitting the Q-distribution, as all three methods are very close in their 
assumptions and limitations) is that it is much less sensitive to the non-flow con- 
tribution compared to the standard method [15]. The sensitivity of these methods 
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to flow fluctuations is not that clear from the above equations. As it follows from 
the simulation results presented in Table 1, the actual sensitivity to fluctuations is 
very similar to that of four-particle cumulant approach. 

Table 1 present the results obtained by different methods on simulated data 
without flow fluctuations, 100% events with v = 0.04 or v = 0.06, and 50/50 
mixture of the above. Last three rows present results from simulations including 
large non-flow correlations (simulated by pairs of particles with erectly the same 
azimuth). The following observations follow: v{2} results, as expected are sensitive 
to both non-flow effects and flow fluctuations (the presented results coincide within 
statistical errors with analytical expectations). «2{4}, «2{BT}, and V2{q-dist} are 
not affected by non-flow. For the case of 50/50 mixture of events with different 
flow values, v 2 {4:} yields value below 0.05 due to flow fluctuations. What was 
somewhat unexpected is that the sensitivity (bias) of V2{BT}, and v 2 {q-dist} to 
flow fluctuations is very similar to that of «2{4}. More detailed examination of 
the Eq. 3 and expansion of the Bessel function around joi up to the second order 
confirms this observation made on simulated data. This is also illustrated in Fig. 2 
(read caption). 

3. Eccentricity and eccentricity fluctuations 

If the initial conditions are totally isotropic - no anisotropic flow is possible. So 
one concludes that at small initial anisotropies the final anisotropic flow must be 
linearly proportional to the initial spatial anisotropy of the system. Although it 
is not clear what one should use (what is most relevant) for the higher harmonic 
flow, in case of the elliptic flow it is relatively straight-forward: one should come 
up with some measure that differentiate the system geometry along the impact 
parameter, in-plane, and, correspondingly, out-of-plane, such as (R y —R x )/(R y +R x ) 
or (Ry — R%)/(Ry + R x ). where radii R x and R y are some characteristics of the 
system in- and out-of- the reaction plane directions. Note, that different definitions 



Table 1. Flow values (in %) obtained by different methods. 





v 2 {2} 


«2{4} 


MBT} 


%{q-dist} 


without non-flow 


100% v = 6.0 
100% v = 4.0 
50/50 v = 4 & v = 6 


6.005±0.003 
3.995±0.004 
5.100±0.004 


6.006±0.006 
3.977±0.002 
4.895±0.007 


6.018±0.004 
4.013±0.008 
4.908±0.004 


5.999±0.002 
3.997±0.003 
4.920±0.004 


including non-flow 


100% v=6.0 
100% v=4.0 
50/50 v=4 & v=6 


6.495±0.003 
4.703±0.004 
5.670±0.005 


6.005±0.007 
4.013±0.002 
4.905±0.008 


6.028±0.004 
4.011±0.002 
4.913±0.005 


6.033±0.004 
4.038±0.009 
4.933±0.005 



G 



S. Voloshin 





1 1 1 , V 1 1 1 1 1 1 1 1 1 


■-. 


i i i i | i i i i . 












BesselJ0(6.*x) 








BesselJ0(4.*x) 








BesselJ0(5.*x) 








0.5*(BesselJO(6.*x)+Besse 


J0(4.K» 






, . 1 , , , . 1 , , . . 1 . 




. , . , 1 . , , , - 



0.35 0.4 0.45 0.5 0.55 0.6 0.65 



Fig. 2. The most left and most right curves shows the expectations for Bessel 
Transform of q— distribution for 4% and 6% flow. The solid black line shows the 
average of these two. The green line shows the expectations for 5% flow. The 
difference in zeroes of black and green lines indicates the effect of flow fluctuations. 



of the eccentricity give different values, but this is not a problem if one uses the 
same definition for all systems/centralities/experiments used in a comparison, and 
as long as the eccentricity itself is small. Recall that the goal of this exercise is to 
exploit the fact that at small eccentricities the final flow should be proportional to 
eccentricity. (Note that such proportionality would not at all be trivial if eccentricity 
is large). Having said all that, the most often used definition of eccentricity is 

{ f~ X 2 , (10) 

where the x direction is taken along the direction of the impact parameter, and the 
average, with some weight, is taken over all points of the system. For the weight one 
can use the participating nucleon or produced particle densities, entropy or energy 
densities, density of binary collisions, etc.. A priory it is not known what weight 
would be the best to use, but again, note that at small eccentricities it should not 
matter provided the same weight is used for all the data. 

In the standard optical Glauber definition of the eccentricity the weight func- 
tion is assumed to be the same for all events with the same impact parameter; any 
event- by-event fluctuations are excluded. In reality, even at fixed impact parameter, 
the eccentricity could fluctuate for many different reasons: the number of individ- 
ual nucleons (quarks, partons) participating in the collision could fluctuate, their 
position in the transverse plane also fluctuate, the multiplicity of produced particles 
in each of "individual collision" could fluctuate along with deposited energy, etc.. 
The first serious attempt to calculate the role of such kind of fluctuations for flow 
analysis in Monte- Carlo Glauber model was done by Snellings and Miller [13] who 
calculate the fluctuation in 

£ std= {y2+x2 y (11) 
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with x, y being the position of the participating nucleons. It was also soon realized, 
that to take the fluctuations fully into account one should take into account also 
fluctuations in the direction of the major axes of the overlap region as well as the 
fluctuation in its center of gravity [17]: 



-part 



(y' 2 



r'2\ 



(12) 



where the eccentricity is calculated relative to the new coordinate system defined 
by the major axis of the initial system region. The results of the corresponding 
calculations are presented in Fig. 3. Note significantly larger difference between the 
standard and participant eccentricities in Cu+Cu collisions compared to Au+Au 
collisions. Using e part instead of e st d significantly improves scaling of i^/e [17]. 

With elliptic flow assumed to be proportional to the eccentricity in an event, 
vi oc £, flow fluctuations should be determined (at least partially) by the fluctuation 
in eccentricity. Thus, «2{2} oc e{2} and t>2{4} oc e{4}, where 

£ 2 }; ( £ {4}) 4 e2( £ {2}) 4 -( £ 4 ). 



(e{2}Y 



(13) 



The results for e{2} and e{4} are also presented in Fig. 3. Two remarks are in 
order for this plot. The ratio of e{2}'s for Cu+Cu and Au+Au collisions is larger 
than unity, which, is used, would improve the scaling observed in [17]. Very small 
values of e{4} for mid-central and central collisions could complicate the analysis of 
«2{4}/e{4} which could become very sensitive to the exact scale of flow fluctuations. 
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Fig. 3. Eccentricity values for Au+Au and Cu+Cu collisions at ^/sjvjv = 
200 GeV as function of impact parameter. For comparison, the solid line in 
the left panel shows the results of optical Glauber calculations used in [12, 10]. 



4. Conclusions 



A new method of flow analysis, the Bessel transform of g-distribution has been 
introduced. It is shown that this method is very similar to the Lee- Yang Zeroes 
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(using sum generating function) and the fitting of ^-distribution methods. It was 
also shown that its sensitivity to flow fluctuations (as well as that of the Lee- Yang 
Zeroes method) is similar to the four particle cumulant results. 

Eccentricity fluctuations relevant for flow measurements v{2} and v{4} are 
presented for Au+Au and Cu+Cu collisions at RHIC energies. 
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